Controlled Manipulation of Individual Vortices in a Superconductor 
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We report controlled local manipulation of single vortices by low temperature magnetic force 
microscope (MFM) in a thin film of superconducting Nb. We are able to position the vortices in 
arbitrary configurations and to measure the distribution of local depinning forces. This technique 
opens up new possibilities for the characterization and use of vortices in superconductors. 



Quantized magnetic flux tubes, called vortices, allow 
superconductivity to survive in high applied magnetic 
field in type-II superconductors, making them more tech- 
nologically relevant than type-L Each vortex has a non- 
superconducting core with a radius on the scale of the 
coherence length £, and a circulating super current that 
generates one quantum of magnetic flux, 3>o = h/2e, 
over the scale of the London penetration depth A. Al- 
though the paired electrons in a superconductor carry 
charge without resistance, a current will exert a mag- 
nus force on all vortices, which results in dissipation, 
if any of them move. Vortex motion is both a chal- 
lenge and an opportunity. The challenge is to under- 
stand and reduce uncontrolled vortex motion. A vortex 
may be pinned in place by co-locating its energetically 
costly non-superconducting core with a defect that lo- 
cally suppresses superconductivity. Decades of materials 
research have characterized pinning strengths and engi- 
neered defects to increase pinning Continued re- 
duction in uncontrolled vortex motion will open up new 
applications, both for quiet circuits in sensing and com- 
munication, and for large currents in high-field magnets 
and power distribution. 

Controlled vortex motion, on the other hand, has great 
prospects for logic applications and for fundamental sci- 
ence. Collectively controlled vortex motion can serve 
as a rectifier, 5 a vortex ratchet mechanism can perform 
clocked logic, 6 and vortices can control spins in an adja- 
cent diluted magnetic semiconductor, 7 while vortices ad- 
jacent to an electron gas in a quantum-Hall state may 
allow the creation of exotic quantum states. 8 A pro- 
posed test£ii£ of the long-standing idea that vortices may 
entangle like polymers^ requires controlled local ma- 
nipulation of single vortices. Vortices are of theoreti- 
cal interest for their own sake j 12 i 13 as clues to the un- 
derlying superconductivity ^ 14 ! 15 as analogs for interact- 
ing bosons^ and as model systems for soft condensed 
matter^ 

Previous experimental manipulations of single vor- 
tices have applied relatively delocalized force a 17 i 18 i 19 i 2Q i 21 
or have not controlled the vortex motion i 22 i 23 Here we 
demonstrate vortex manipulation with nanoscale control 
and show that we can quantitatively measure the local 
depinning force. In magnetic force microscopy (MFM), 
the sample exerts a measurable force on a cantilever with 
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FIG. 1: Manipulation of vortices to spell SU. Colorbars give 
Af [left for (a-c), right for (d)]. (a) An unmanipulated con- 
figuration of pinned vortices after initial cooling to T — 7. OK, 
imaged at a scan height z = 300nm. (b-c) Intermediate con- 
figurations after manipulation of some vortices in the temper- 
ature range 7.0 to 7.2K, imaged at z = 300nm. (d) Final con- 
figuration after completing the vortex manipulation at 7.2K, 
imaged at z = 120nm and T = 5.5K for better resolution and 
stronger pinning. 



a sharp magnetic tip, such that scanning the cantilever 
at a constant height z above the sample provides a map 
of magnetic features. MFM has been used for a variety 
of vortex experiments j 22 i 23 i 24 i 25 i 26 i 27 Many experiments 
image the cantilever deflection, which is proportional to 
the vertical force, F z . For improved signal-to-noise ratio, 
we use frequency modulation (FM) mode, 28 in which the 
imaging signal is a shift Af in the cantilever's resonant 
frequency, /q. The images show the variations in the 
derivative of F z , dF z /dz = — 2fcA///o, where k is the 
cantilever's spring constant* 2 ^ We use the lateral com- 
ponent of the force Fi at = \F x x + F y y\ to pull or push 
vortices. 

We used a home-built variable-temperature MFM^S to 
study a 300nm thick Nb film sputtered onto a silicon 
substrate*^ The midpoint transition temperature is T c = 
8.6K with a transition width AT C = 0.6K, as measured 
by magnetic susceptibility, and A = 90nm. 32 

Figure [1] shows controlled vortex manipulation. We 
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cooled the sample to T = 7K in an external field of a 
few Gauss with a polarity giving an attractive tip- vortex 
force. Figure shows a disordered arrangement of vor- 
tices pinned in the sample. The lack of observed vor- 
tex motion during the scan indicates that at this height 
of z = 300nm we are imaging vortices without depin- 
ning them. To move a vortex, we position the cantilever 
tip over it, reduce the oscillation amplitude, descend 
to z ~ lOnm, move parallel to the surface, and then 
withdraw vertically, leaving the vortex at a new loca- 
tion. Some locations require several attempts at different 
heights or lateral offsets to move the vortex, or different 
pulling directions and slightly elevated temperature to 
reduce the pinning force. The exact final location of the 
vortex must depend on the local pinning potential; pin- 
ning sites appear to be dense in the Nb film. Images 
(Fig.[T]>d), each taken after several manipulations, show 
the successful repositioning of vortices to write SU. This 
procedure is akin to atom manipulation in scanning tun- 
neling microscopy^ 

We also quantified the forces required to depin vortices 
at 5.5K. Figure [2] shows a series of images, with decreas- 
ing z, of 8 isolated vortices at a vortex density corre- 
sponding to 4G with a polarity giving repulsive tip- vortex 
force. As z decreases, the increased dF z /dz leads to in- 
creased signal strength, while the increasing Fi at depins 
vortices. Vortices at different locations depin at differ- 
ent heights, indicating a distribution of depinning forces. 
In this case, F\ at is weak enough that a depinned vortex 
simply finds a better pinning site nearby, but we have 
also applied larger forces to sweep the field of view clean 
of vortices. 
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FIG. 2: Vortices in a Nb film imaged at T — 5.5K as a function 
of decreasing scan height z. (a-e) Images at 5 different values 
of z. Red right arrows indicate vortices depinning for the first 
time. A single vortex remains unmoved even at the largest 
applied force in (e) (yellow left arrow), (f-i) Cross-sections 
along the crosshairs in (a) for the first four values. The lines 
are from the monopole-monopole model, (j) Closeup on the 
boxed region in (e). 



We fit the data in Fig. as well as a second dataset of 
31 vortices at a vortex density corresponding to 17G, also 
at 5.5K. We model the cantilever tip and each vortex as 
monopoles* 2 ^ Cuts through the fits are shown in Fig. [Sff- 
i). The resulting max(i^) (Fig. [3^) has systematic error 
bars of less than 50% from the modeling and from the 
uncertainty in k. These systematic errors do not affect 
the relative distribution of depinning forces. 

Most vortices clearly moved by visual inspection dur- 
ing one or more scans. To quantify our detection thresh- 
old, we bootstra p 29 ! 34 the entire 8-vortex data set and 
one quadrant of the 31- vortex data set.— The upper limit 
to undetected vortex motion at the 95% confidence level, 
5R, is shown in Fig. [3)3. The stationary vortex in Fig. [2] 
did not move by more than 9nm. 

Identifying the z-value at which each vortex moved 
gives a histogram of vertical depinning forces (Fig. [3b) 
ranging from 12pN to more than 32pN. Since it is ac- 
tually Fi a t that causes the depinning, we use modeling 
to estimate it based on the measured max(i^). Depend- 
ing on the model, the ratio between the maximum F\ at 
and the maximum F z ranges from 0.3 for a pyramidal tip 
less sharp than ours 35 to 2/3^ « 0.38 in the monopole- 
monopole model. Using 0.38, the observed lateral depin- 
ning force ranges from 4pN to above 12pN, or 15pN//im 
to 40pN//im normalized by film thickness. This tech- 
nique could be advanced further by using vertical can- 
tilevers with lower spring constants to directly measure 
the lateral force. 

The manufacturer-specified critical current for similar 
films is 50 ± 30mA per transverse /am at 4.2K, 31 equiva- 
lent to a depinning force of 104±69pN. Other single vor- 
tex pinning measurements in Nb used a transport current 
to supply a forc o 17 i 18 i 19 i 2Q and various stationary probes 
to detect motion. For comparison, we normalize previous 
results by the film thickness and extrapolate to our re- 
duced temperature using the power law F ex (1 — T/T c ) 7 , 
where 7 is an experimentally determined exponent that 
varies greatly between experiments. The inferred depin- 
ning forces are 123pN//im, 20 80pN//miji£ 58pN//imji£ 
5pN//imji£ and 0.6pN//im^ These experiments have 
motion detection thresholds ranging from a few hundred 
nanometers to microns, except for [19[, which reports a 
resolution of ~ 20 — 40nm, achieved with an array of 
stationary Hall probes. 

Our best threshold for vortex motion detection is bet- 
ter than lOnm and is limited by the signal-to-noise ratio. 
The characteristic scale for changes in the pinning po- 
tential is the coherence length £ w 10 — 20nm j 17 i 25 We 
can therefore detect all vortex depinning events with a 
quantitative determination of the locally applied depin- 
ning force. Imaging the vortices before, during, and after 
the depinning has great prospects for correlating pinning 
with topography, for determining the pinning landscape 
directly, and for studying single- vortex dynamics. 
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FIG. 3: (a) Maximum of the vertical force applied to a vortex, max(F z ), as a function of the measured tip-sample separation, 
z, and the modeled monopole-monopole separation, z + A + d, for two data sets taken with the same cantilever. Solid green 
lines indicate the systematic error associated with modeling the data. Dashed purple lines also include the uncertainty from 
the cantilever spring constant, (b) The threshold for detecting vortex motion at the 95% confidence level. The 4G dataset 
has better resolution because the data was less noisy. The curves are from the monopole-monopole model.— (c) Histogram 
of number of vortices depinned vs. max(i ? 2 ). The width of the histogram bars is the difference in force between heights of 
successive scans. The arrow on the last blue bin indicates that a single vortex did not move at the largest applied force, as 
shown in Fig. [2] 
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Supplemental Material 



I. MFM TIP FOR FIG. CD 



We acquired the data shown in Fig. [T] using a MikroMasch silicon cantilever with fo = 81kHz, k = 4.0 ± 0.6N/m, 
and an electron beam deposited (EBD) tip, with an iron film of nominal thickness 40nm evaporated onto one side of 
the tip. 



II. MFM TIP FOR FIGS. f2l f3l ISll 



We acquired the data shown in Fig. [21 and analyzed in Figs. [3] and [HU using a Veeco Instruments silicon cantilever 
(model MESP) with /o = 71kHz, k = 2.1 ± 0.7N/m, and a nominally uniform CoCr Veeco-proprietary coating. 



III. ANALYSIS 



The field lines for an isolated vortex in a film of finite thickness t and infinite width are given in For z,t ^> A, the 
vortex field lines are those of a monopole of strength 2$o//io located at Ro = (#o, Vo) a distance A below the surface, 2 
where /io is the magnetic permeability of vaccum. We model the tip as a monopole m offset from the physical tip of 
the pyramid by a distance d, 3,4 as shown in the inset to Fig. ISTk . This monopole - monopole model gives a frequency 
shift: 



fo „ d£J ortex f m$ —(R ~ Ro) 2 + 2(z + A + d) 2 



2k dz 2k 2tt 



5/2 

(R-R ) 2 + (z + \ + dY 



(SI) 



per vortex, where £?j ortex is the vertical component of the magnetic field from a vortex and R = (x,y). We fit an 
image to a sum of terms A/ m , one for each vortex, plus d/offset, a z-dependent, ^-independent, offset due primarily 
to the laterally uniform Meissner repulsion of the tip from the sample. The free parameters of the fit are Ro for each 
vortex and df offset, font /2k, d + A, all the same for all vortices in the image. The latter two fit parameters (Fig. ISTk.b) 
varied by about 25% over the measured z range, perhaps because the signal is increasingly dominated by the very tip 
of the pyramid for smaller z. The fit value of fofh/2k corresponds to a monopole moment on the tip equivalent to 
about 5 vortices. 

We also tried a dipole model for the tip, and found that the monopole model gave smaller x 2 an d more consistent 
values for the tip parameters. Possible refinements include modeling the tip as an extended object, or modeling 
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FIG. SI: (a) Fit parameter d + A from the monopole- monopole model, averaged over all as yet unmoved vortices at a given 
height, vs. z. Inset: sketch motivating the monopole model for the tip. (b) Fit parameter forh/2k, also averaged over all as 
yet unmoved vortices at a given height, vs. z. (c) The maximum vertical force gradient majt(dF z /dz) at each height, plotted 
against the measured z and against z + A + d, with A + d obtained from the fit at the lowest z- value. The latter plot shows the 
expected power law behavior with p = 3; the p- value for each pair of points is indicated. 
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the vortex using the full expressions for a film of finite thickness. However, both the peak frequency shift, which 
gives the maximum vertical force gradient, and the vortex center, which enables us to identify vortex motion events, 
are already well determined by the monopole-monopole model. Also, the bootstrapping required to determine the 
confidence intervals involves hundreds of iterations of simultaneous multi- vortex fits, and would become cumbersome 
with more complicated models. 

The peak frequency shift, averaged over all as yet unmoved vortices at each height, gives the maximum vertical 
force gradients max(<9i^/<9z) (Fig. ISTb). If the monopole-monopole model is strictly valid with no z dependence of 
the parameters, we expect max(dF z /dz) oc (z J r\ J rd)~ p with p = 3. The value of p for fits to each pair of consecutive 
ma,x(dF z /dz) points ranges from 2.95 to 3.03 (Fig. ISTb), indicating good agreement. The measured max(<9F 2 /<9z) 
data, extrapolated to infinity with p — 2.95 and integrated from infinity to z, gives the maximum vertical force 
max(F 2 ) (Fig. 33*). 

To characterize the possible systematic error associated with this extrapolation, we note that max(dF z /dz) plotted 
vs. zona log-log plot has negative curvature. We extrapolate using ma,x(d F z /dz) oc z~ xm for an upper bound, and 
ma,x(dF z /dz) oc z~ 4 , the power expected for a dipole tip, as a lower bound. This systematic error could be reduced in 
future experiments by taking additional data at larger z. The other main source of systematic error is the uncertainty 
in the cantilever spring constant. 



A. Detection threshold for vortex motion 



The threshold for vortex motion detection is determined by bootstrapping an 8- vortex dataset at 4G and a 9- vortex 
dataset at 17G. In each case, we fit all 8 or 9 vortices simultaneously to Eq. [HH allowing forh/2k, d + A, xq and 
2/o to vary separately for each vortex. Because fork /2k and d + A are properties of the tip, we expect them to be 
constant for all vortices at a given height; indeed we find that fom/2k varies by less than 2%, while d + A varies 
by less than 5%. We use at least 600 bootstrapping iterations to find the statistical distribution of (xo, yo) pairs 
for each vortex center at each height. Each vortex has a bootstrapped set of center parameters XQ nvi and yo n v i 
where n labels the scan height, v labels the particular vortex, and the bootstrap iteration index i runs from 1 to at 
least 600. To find the confidence levels for the motion detection threshold (Fig. O)), we looked at the distribution of 

SR = \l ( X 0n,v,i ~ ^0 n +l,^,i) 2 + (yo n ,v,i ~ ^0 n+ l, v ,i) 2 . 

We have also derived a formula for the detection threshold. For concreteness let us assume that the tip-vortex 
interaction is indeed described by the monopole-monopole model. The motion detection threshold is proportional to 
the ratio between noise of the measurement (SAf) and the slope of the measured curve: 

SR - 2k wmmm - (S21 

Assuming that the vortex first moves when the lateral force is maximal, which in the monopole-monopole model 
occurs for R* = (z + A + d)/y/2, we find: 

SR = A(z + \ + d) 4 , (S3) 

where A = fp=^ (SAf/fo). Equation [S3] is the formula plotted in Fig. [3jo, where we used A + d from a linear fit to 
Fig. 15Tb and proportionality constants, A, from fitting the data in Fig. [3)3 to Eq. fS3l 
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